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and independent of depth 
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in-phase component of the tidal height 
out-of-phase component of tidal height 


forcing function 


1, Introduction. 

The objectives of this study are to develop a numerical 
method for estimating tidal constituents in the open ocean using 
fundamental theoretical concepts, to adapt the method to com- 
puter use, and to apply it to the analysis of the tidal regime 
in the Gulf of Mexico, 

The ability to predict tidal behavior in the open sea is 
useful in oceanographic research and coastal engineering. It 
permits removal of tidal effects from deep-sea observations 
and prediction of the effect upon the tides of alterations of 
the coast line. 

Present tidal knowledge is largely limited to the results 
of harmonic analysis at selected reference stations and empirical 
Jag and amplitude corrections for closely adjacent stations. This 
technique is, of course, limited to the coast, but the work of 
Dietrich [1] , subjectively fitting a two-constituent (M2 & K1) 
model (see below) to island and coastal date between coastlines 
extends the results into the open ocean. The accuracy of 
Dietrich's work is considered good in dense data regions, such 
as the North Atlantic; but in sparse data areas no general 
agreement has been reached concerning tidal regimes. The dis- 
tortion of tides in shallow water cannot be accurately assessed. 
Thus the results of the method are, to some degree, in error. 

A more satisfying approach is to apply the fundamental 


hydrodynamic equations and calculated luni-solar forcing function 
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to as realistic an ocean basin as possible. 

Two methods of attack are the integration of the basic 
equations and the approximation of the equations by finite dif- 
ferences, leading to a solution by numerical methods. In either 
case, a boundary value problem results in which increased real. 
ism is accompanied by increased mathematical difficulty, where- 
as oversimplification can render the problem physically meaning- 
less. In this paper discussion will be confined to numerical 
approaches. Hansen B used finite differences in his treatment 
of the M2 tidal constituent in the North Sea with excellent 
results, both in tidal heights and currents, as shown in Figure 
1. Friction was included and the beta-plane approximation was 
used. Rossiter's approach [3l involves the use of spherical 
coordinates and frictionless flow in a meridionally bounded 
ocean of constant depth. Coriolis effects necessitate spec- 
ial care near the poles and the parameters were modified to 
insure convergence in these areas. In another application lat- 
itudinal boundaries were imposed to approximate the Atlantic 
Ocean. Relaxation techniques were employed successfully by 
both Hansen and Rossiter. Their results and Rossiter's recom- 
mendation that the method be applied to irregular shaped basins 


were important motivations for the present study. 
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2. Tides in the Gulf of Mexico 

The tides in the Gulf of Mexico are modest in range, but 
diverse in character. Along the Gulf Coast the tides change 
abruptiy and assume many forms. Figure 2 illustrates the 
distribution of these regimes. Along the entire Atlantic coast 
of the United States tides vary widely in range, but there is 
no variation in semi-diurnal character. 

In the light cf this constancy the abrupt changes in the 
Gulf are very dramatic. One example of considerable variation 
of the tidal regime is in the northeast portion of the Gulf. 
At Cedar Keys, ἃ mixed, predominantly semi-diurnal regime 
exists. The tidal curves vary markedly in both directions along 
the coast. At Pensacola, 200 miles west of Cedar Keys, the 
diurnal influences dominate. Similar changes are encountered 
in other areas along the Gulf Coast. A comprehensive descrip- 
tion of the tides along the Gulf Coast and the criteria for 
regime classification have been given by Marmer [1] . Many 
explanations have been advanced for these shifts in the basic 
form of the tides. Harris postulated in 1897 that the semi- 
diurnal tides were the extension of the Atlantic regime, while 
the diurnal tides resulted from the natural period of the Gulf 
basin. 

Grace [5] compiled the various theories and applied an 
inventive analytical technique in 1932. He divided the Gulf into 


five sections in each of which the depth was treated as a 


TT 


constant. The linearized, inviscid hydrodynamic equations of 
motion were solved within each section. Continuity was satis- 
fied a the interior boundaries between sections and the bound- 
ary conditions of no normal flux were met at the coasts. The 
studies resulted in a unique set of corange and cotidal lines. 
Based on this work Grace put forth a theory of tidal dynamics 
in the Gulf: the semidiurnal tidal wave enters the Straits of 
Florida, moves around the Gulf coast in a contra solem amphi- 
dromy and exits through the Yucatan Strait six hours later; 
the diurnal tidal wave enters the Straits of Florida and prop- 
agates in such a way as to arrive at most points along the 
coast nearly simultaneously. 

Several topographical features of the Gulf may well have 
discernable effects on the tides. The continental slope off 
the west coast of Florida is very steep and could act as a 
reflective surface to an incoming tidal wave. The intrusion 
of the Desoto Canyon toward the northern coast and the 
Mexican Basin are both relative deeps with unknown effects on 
the natural period of the Gulf. The Florida Keys and the 
broad shelf off the west coast of Florida could mask influences 


of tidal waves entering the Straits of Florida. 
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3. The Theoretical Development 


Introduction 
To describe the tides utilizing the basic hydrodynamic 


equations one must resort to describing the separate tidal 
constituents. This limits each solution to a particular known 
period and provides for a complete solution through a summation 
of constituents. Except in specialized cases this summation 
limits the problem to linearized equations of motion. The work 
described below is directed towards a solution for the principal 
lunar semidiurnal constituent, Me. Coastal observations indicate 
that knowledge of this constituent should permit explanation of 
these observations and description of regimes in the open sea. 
Harmonic analyses presented by Marmer indicate that M2 and 
K1, the lunisolar diurnal constituent, are the dominant constit- 
uents. The variability of the amplitude of K1 is from .59 

feet to .09 feet, ignoring the results of short-term observa- 
tions. ‘The amplitude of Me varies from zero to 1.07 feet by 
the same criterion. Thus the spacial distribution of the M2 
constituent should be indicative of the distribution of the forms 


of the tides. 
The Equation of Motion 


From the equations of motion and continuity we may derive 
a second-order elliptic partial differential equation relating tidal 
height, depth, and their derivatives. The analysis is performed 
on a beta-plane. Advective and frictional terms in the equation 


of motion are neglected. The hydrostatic approximation is made. 
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The beta-plane assumption restricts the size of the 
basin in which the results will apply, but is considered accep- 


table in the Gulf, 


Friction is present and does influence tidal distributions, 
especially in shallow water. Frictional effects are, however, 
believed to be highly non-linear, thus further complicating the 
equations. (The work of Dronkers| 6] in linearizing this term 
for shallow water is an example of these complications.) Fur- 
ther justification for neglecting friction is that only the per- 
imeter is in shallow water while most of the interior is in 
relatively deep water. 


With the stated assumptions the dynamic equations are: 


E PAEA 
ge -FV= 98% 


ON + Fu gi (1) 


where, 


u is the velocity component in the x direction (positive 


East), 


v is the velocity component in the y direction (positive 


North), 
F is the Coriolis parameter, 
g is the acceleration of gravity, 


f is the tidal height above undisturbed sea level (which 
is mean low water) and is positive upwards, and 
J^ is the forcing function (see below). 


The hydrostatic approximation implies that u and v are 


independent of depth. 
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The equation of continuity 16: 


E NV), Om 
Эмо). 9f v) <š =o 


where h is the basin depth. 


(2) 


Since the frequency of the desired harmonic constituent 
is known, sinusoidal forms for the dependent variables are 


assumed: 


fers Reli Y) ae 

U = Ref u(x, ner (3) 
М = Ref v (x, vY e? й, 

£= H cos OT 


where the notation, ке) $ indicates the real part of the 
quantity in brackets, 

O” is the angular frequency of the constituent under exam- 
ination, and 

H is the amplitude of the forcing function. 

Solving the equations ef motion for u and v and substitut- 
ing the results into the continuity equation yields one partial 
differential equation for the tidal height. Writing J = ge 4 72. 
(where, for е is the component in phase with the 
forcing function) and separating the differential equation into 


real and imaginary parts gives 


(AA UON S. , Oh OH 
^ (41-4) * FR SR - & BY SE "8484 ) 
= Së Gr TË K = o 
Bs аа 3 па) au (4) 
on O Š x E eh CH eu oJ 
hv Ji *3x Ox. .' G O4 Ox “9424 
едм (н) СЕ = О 
T Өх Ta 


х 
-ᾱ- 
с 
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where 7° = д^ の < 


2x*・ Sut 
Rossiter used spherical coordinates in treating basins 
covering large areas of the earth. Thus his equations differ 
in detail at this point although the form is very similar. 


Hansen's work in the North Sea included friction. The 


added terms in the equations of motion had the form, 


Fe: KU and ky: KV (5) 


where R is a suitably chosen proportionality factor. The 
excellence of his results is a strong argument for this some- 
what ad hoc hypothesis. Tn the situation involving considerably 
greater bottom relief this form is inadequate due to the var- 
iability of the friction effect with depth. In deep water 
friction will not influence tidal distributions as much as bottom 
topography which enters the equations through continuity con- 
siderations. 
Boundary Conditions 

The boundary conditions determine the velocities normal 
to the boundaries. At solid boundaries these velocities must 
be zero. At the passages into the Gulf values must be assign- 
ed to the associated tidal currents. Following the treatment of 


tidal heights, we write 


U = UK| エキ 上 ん いい 2 
(6) 


V= VKI i VET 
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Substituting these relations into the solutions for u and v from 
the equations of motion and separating into real and imaginary 


parts gives 


VKI = - ges е” 


a ә(Ж-ндеса 2 
Уқ2- 26<(6%-ғ?) jo 20 ) ў 
Se ο 
ur zzi]? 25 pa ع‎ 
D-H). 25 
UK2= zr gl i FS. 


Observations of tidal currents to determine their con- 


ja e Об») 


stituents are not frequently made. For accuracy continuous 
observations are required over long periods of time. Tidal 
current measurements made by Harris in 1887 and supported 
by Grace's work were used in the Yucatan and Florida Straits 
in this study. The accuracy of these calculations is question- 
able but the order of magnitude is correct. By allowing flow 
through these channels the effects of a co-oscillating tide are 
included in the solution. 
The Forcing Function 

The forcing function is basically the gravitational effect 
of the ideal body associated with the constituent to be examin- 
ed. It is obtained from equilibrium tidal theory (for example 


see Proudman[7]) in the form, 


të cos? P eos“ d Cos を ke 





る 
Mm 
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for the lunar semidiurnal constituents, and 


A Mn (3 е 2% сос 2. d 685 LE 
H=ald) me Е и 
for the lunar diurnal constituents. 
Flere, 
a is the radius of the Barth 
d is the distance from the Earth to the Moon 
Mm is the mass of the Moon 
Me is the mass of the Barth 
Ф is the latitude of the point under consideration 
di is the declination of the ideal body 
Y is the longitude relative to the reference longitude. 

For computations the central meridian of the Gulf, 90°w, 
will be used as a reference. 

Similar expressions apply for solar diurnal and semidiurnal 
constituents and also for those with periods of up to two weeks. 
In this work these constituents have relatively small amplitudes 
and so are not treated. 

Over the area of the Gulf the lunar semidiurnal forcing 


function varies by about eighteen percent. 
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ho Numerical "Techniques 
Finite Difference Approximations 
To put equations (4) into & form more suitable for num- 


erical treatment the derivatives are approximated by central 


differences, 
Oz ー Ἔχει - Ἔχ-ι の 4 - Zw - Ëv-i 
д Х /о LE ; OY O Ze (10) 


and 


| 
(v'z), A: ἐν + Eye Fey 2 - | (11) 


whereZ is a dummy variable and the subscripts refer to loca- 


tion in the grid arrangement as shown below,. 
Yai 


[ 


= CE 


Kl 5 


K+, 


Y-1 
Central differences utilize values sampled at an interval of two 


grid distances. This has a smoothing effect on the data. 
Because the bathymetry is normally given to the nearest fathom 
smoothing is beneficial in decreasing errors in depth gradient 
se to sounding inaccuracies. The tidal heights are much smaller 
than the depths. Because of this, solutions are very sensitive to 
sounding errors. For this reason central differences are used 
in all finite difference енен 

The use of central differences at a point on the bound- 
ary introduces välues for tidal heights outside the boundary into 
the equations. 

For an interior point, equations (A4) with finite differences 


become 
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K^ 
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ь н е (пеку) + (3), [e 2007 = o. 
Consideration was given to variable grid distancing’ to permit 
detailed analysis in areas of large depth gradients. However, 
the stability of the computation about a point with different 
distances to each of the adjacent points is questionable, and 
the finite difference form of the equations is quite complex. 
Appendix B contains the finite difference equations for a 
variable grid distance. The density of bathymetric data in 


the Gulf is such that a grid finer than that used below would 
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require considerable interpolation and thus no improvement in 
the solution would appear. 

Rossiter used a finite difference equation involving der- 
ivatives at the boundary and three interior points to produce 
the tidal heights outside the boundary. Although it is valid, 
Rossiter's procedure has the disadvantage of introducing para- 
meters associated with non-adjacent points into the equations. 
This complicates the system of equations in a relaxation tech- 
nique by increasing the number of points involved in each equa- 
tion by two. It also disrupts the iterative form designed for 
computer use, which wil be discussed later. 

At a boundary point, the basic equations (12), derived 
from continuity considerations and the applicable boundary con- 
ditions are both in effect. These boundary conditions manifest 
themselves in the forms of prescribed velocities normal to the 
boundaries. 


The velocity equations (7)in finite difference form are 


VK | = ey r (る ん ыі: E. (7-4) ў 
Мк у (ee -( ση, ο. (Aha) | 


(13) 


əкі A Ге ею) 
оке (A TC te [CR 


al 


Thus, for a given normal velocity, the tidal height for the 
point outside the boundary may be solved for in terms of three 
interior tidal heights. Substituting these solutions into the 
equation for an interior point now incorporates the dynamic 
equations, the continuity equation and the boundary condition 
into a set of two equations at the point, as well as eliminate 
ing values outside the boundary from the calculations. 

The substitutions and resultant equations will be unique 
for each boundary or corner. As an example, to solve for a 
point on the south boundary with a known velocity, equations 
(13) are used in the form, 

る 
(Aa = (SH) yr GLa dear Ue ЕЗ) 
(14) 
re) + GE) 
γ-ι 2/ С A а he, 

and substituted into the interior equations (12) to eliminate 


values involving the exterior point Y-1. 


This results in 


(A) ya [Bho] + (32) [5 4h0] + (7), (4) + 
+ (Lx L4ho + Che he - E) | 
+ (FH) x LAW (serine (14707) ]+ 4h. H 15) 
+ ア ). [ 569. leho] Е ve konano 
T мка Со (Δ 11ο - hasce haa) 
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ес 


AR [She] と ( TA) s [E Aho | ан) (4 ehe) 
D A 4h。 тача = д") | 

+ (37) pal Ano Сака па (А ZEN 

СЯ. кч = 

τω. E (a No” Nur н ү 
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To accommodate the boundary conditions involved at a corner 
values of u and v must enter the basic equations (12). This 
is accomplished by simultaneous solution of equations (13) to 
eliminate terms involving the two heights at the same point 


from the combined equations. The results are 


E QE pe ока [2658 - VK ене 

(Кн) а ы UK| Ec анеки RE 
κε と と < ) 

(Ж + а (σι) K-1 + VK2 eel UKI Er 


(16) 


GH ue = (FH з), с KESH: + әс Zei) 


Combinations of these equations will meet any corner boundary 


condition. For example, substituting to eliminate values at the 
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point X+1 in the south boundary equations will give equations 
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The forms of the equations for each boundary and corner 





point are presented in Appendix A. 
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Methods of Solution 

When equations (12) and their equivalents incorporating 
boundary conditions (see Appendix A) are applied at each point 
of a grid consisting of n points describing the body of water 
studied, à system of 2n equations in Zn unknown tidal heights 
results. Both methods of solution undertaken here are pred- 
icated upon the presence of a dominant diagonal in the system's 
coefficient matrix. That is, in the application of the equations 
to each point there is one coefficient much larger than all 
others in the equation and in each successive equation the 
largest coefficient multiplies a different unknown value of the 
tidal height. Both methods are designed to obtain a solution 
by altering the unknown tidal height associated with the largest 
coefficient. 

Gauss - Seidel iteration (see for example, Hildebrand [8] ) 
is the method used in the computer program. Briefly, this 
method consists of obtaining new values for the variable with 


the dominant coefficient, say X4 , in an equation of the form, 


A,X, + A2 Xa + Я An Xn = О (2120) 


by the formula 


N 
| 
X, ----2,, Ам Хм, 


This satisfies the equation with approximate values of Хуу... 


(18) 


Ха by varying Xjonly. This procedure is carried out for each 


equation in the system following a first guess of the value of 
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each variable. Repeated use of equation (18) usually brings 
about "convergence", in the sense that each of the equations 
is satisfied by the same set of unknowns to within prescribed 
tolerances. Gauss - Seidel iteration is not highly sensitive 
to the first guess values. 

Relaxation is a similar process but permits selective 
variation of the unknowns to bring about convergence. Using 


the same example equation the first guess will give 


E SW CAES uec E Αν Κω = R 


where R is the "residual", which should vanish to achieve con- 
vergence of the system of equations. Changes will be made 
in the unknown associated with the largest coefficient following 


the rule 


/ 
Lt hte ae - 


where N is the "Relaxation Coefficient" and may take on values 


in the range 
= (20) 


IMA 
When “ =1 the residual wil be reduced to zero and the 
process is very similar to the iteration described above. When 
>| the residual tends to alternate sign on successive iterations 
and decrease in magnitude. This process usually speeds up 
convergence and is called overrelaxation. When!<\the residual 
tends to be reduced in magnitude and keep the same sign thus 


approaching zero from one side rather than oscillating around 


zero. ‘This is called underrelaxation. 
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Su Hand Calculation of the M2 Tidal Constituent in the Gulf 


To gain an understanding of the behavior of a tidal solu- 
tion in the Gulf and to supplement the computer work a hand 
solution was undertaken for a relatively coarse grid. The grid 
arrangement consists of 16 points with a grid distance of 180 
nautical miles. Figure 3 shows the arrangement. This grid 
represents the Gulf fairly wel, but has the disadvantage that, 
in areas of large depth gradient, the terms in the equations 
may not occur with the relative magnitudes desired for a rel- 
axation solution. 

Applying equations (12) and appropriate boundary conditions 
to each of the 16 grid points produces a system of 32 equations 
for 32 unknown tidal heights. Of these, eight do not contain 
the dominant coefficients desired for relaxation. Six of these 
eight are for interior points with eight terms in each equation. 
This means that successful reduction of the residuals of these 
equations, especially in the later stages of the solution, may 
depend upon systematic adjustment of the number of tidal heights 
to avoid deleterious effects on the residuals at adjacent points. 

In every iteration values of the tidal heights used were 
the results of the last calculations made at the point. Some 
systems (e.g., the method of Richardson as described by 
Hildebrand [з ыа.) ) are cyclic in nature using only calculations 


from previous iterations. The approach in this work is not 
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cyclic，but is based upon maximum reduction of the largest 
residuals consistent with the long-range effect upon adjacent 
points. 

Initial passes were made based on a first guess field de- 
rived from Grace's results and observational data reported by 
Marmer. Over-relaxation techniques were utilized with IN tak- 
ing values from 1,5 to 1.2. The chief drawback of this tech- 
nique is that it can produce oscillations in the values of the 
tidal heights and thus fail to converge. At the first sign of 
an undamped oscillation, I“ was lowered to 1. On subsequent 
iterations the relaxation factor was varied to meet the needs 
of each equation. To determine corrections late in the cal- 
culations systems of simultaneous equations were solved for 
the appropriate values of TN. At this stage the equations 
were very sensitive, and large corrections frequently had del- 
eterious effects on adjacent points. To most efficiently deal 
with these situations tabulations of present residuals and the 
effects of proposed changes in tidal heights upon them were 
necessary, These were most helpful in determining the inter- 
action of changes in the tidal heights. After twenty iterations 
component tidal heights ( апа @) at all interior and solid bound- 
ary points were defined to an accuracy of plus or minus three 
centimeters. Those at the Yucatan Passage and Straits of 
Florida were extremely sensitive to values of input current, 


fluctuating approximately ten centimeters with changes of only 
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one centimeter per second in the associated velocity component. 
The effects of these points on other equations is not large. 
Since the current velocities utilized are not considered suffic- 
iently accurate to warrant it, these points were not relaxed 
to the degree applied to the other points in the grid. 

The results of these calculations are contained in Appendix 
D. Іп general, this solution was very similar to that accom- 
plished using the computer on the same i6 point grid. The 
maximum difference in magnitude of tidal height was 9.06 cms. 
The average difference was 3.23 cms. This computer solution 


is discussed in more detail below. 
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6. Solution by Digital Computer 

Adaptation of the finite difference formulae (l2) to the 
computer is a difficult task. For efficient iteration a proced- 
ure is necessary for expressing all equations in a common for- 
mat. This procedure permits the use of indexing operations 
in the identification and location of equations. This utilization 
of indexing is accomplished by the inclusion of identification of 
adjacent points and their relative positions in the input data. 
Through this, and a coded method for points outside the bound- 
ary the program identifies interior and specific boundary loca- 
tions and simultaneously applies the proper equations to them. 
Fortran 60 program language is used. Based on the needs of 
this particular application, convenient computer capacity, and 
the concept of versatility, the maximum number of points in 
the grid is 1200. These points are numbered serially, and 
points outside the boundary are identified by numbers greater 
than 1200. Since the requirement for ficticious tidal heights 
outside the boundary has been eliminated (see above), the only 
data required outside the boundary are tidal current velocities 
and basin depths. 

The presence of an identification number greater than 
1200 then specifies a particular boundary and two such numbers 
specify a particular corner. ‘The current velocities included as 
data for points outside the boundaries replace tidal heights for 


those points in the equations at the boundary point. This 
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requires alteration of the coefficients in the equations to fit 
the location of the point at which the equations are being ap- 


plied. This is accomplished from the general forms, 


las [αι G за -Ὃι (ο)... 
+ 22 ( 37-H)x- +0: (2) х1 


«Ακτή ・b5( Five » 
+ б.а (57-н) хн - pa ( T2) «4 
n + CJ/ Cz 
(12; = Γαι ( 72) και + bi (SH) xi 
+01 (32)x-1 - ら < ( (CS%-H) xa 
каз (У) - Оз (Si -H)yrı 
+ 0.4 ( 3^ «i + ра (н) 1 /cz 、 


Thus the coefficients depend upon the location of the point. 
Appendix A presents these coefficients in detail and the com- 
puter program is given in Appendix C. 

The system of equations derived from the grid utilized 
in the hand calculations was then solved using the computer. 
Actual depths and currents were employed. This system con- 
verged in 37 iterations.. 

In a study of this nature results are best displayed by 
the analysis of the solution to obtain cotidal and corange lines. 
Cotidal lines are lines along which the high tide occurs simul- 
taneously.  Corange lines are lines along which the tidal ampli- 
tudes are equal.As set forth above these results apply to the 
M2 constituent, and the phases are relative to the passage of 


the moon through the 90th west meridian. The lack of detail 


al 


in this solution causes the corange and cotidal analysis to be 
somewhat subjective. The analysis of the 16 point solution is 
illustrated in Figure i which is based upon coastal data as well 
as the computer solution. It is evident that most of the 
distortion of the tides occurs in the coastal regions. This 
does not mean that the solution is in error. The lack of 

detail does not permit identification of a small-scale amphidramy 
such as that proposed by Grace or seen in the solution discussed 
below. The path of the progressive tidal wave indicated by the 
cotidal lines is logical and amplitude distributions agree with 
observed coastal data. 

To obtain more detail in the solution a finer grid was 
desired. Using a grid distance of thirty nautical miles a grid 
was then constructed in an irregular shape closely approximating 
the Gulf. Real depths were used, taken from U. S. Coast 
and Geodetic Survey Chart 1007. The grid consisted of 504 
points. Current data mentioned above were used in the Straits 
of Florida and Yucatan Strait. The resulting system of equa- 
tions failed to converge after 2500 iterations. In the initial 
iterations the tidal heights increased slowly in magnitude. 

After 25 iterations the solutions were unrealistic and they 
increased more rapidly in magnitude. Computational instability 
seemed to be present and further efforts were directed at a 
more tractable situation. 


Using the same grid 8 constant depth of 1000 fathoms 
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was introduced. The system of equations converged in 378 
iterations with a relative convergence criterion of 0.1. That 
is, each tidal height within the grid must change by less than 
10% of that of the previous iteration before the iterative 
procedure is stopped. The maximum change of tidal height 
over the last three iterations was .12 centimeters. These 
results are tabulated in Appendix D. Figure 5 shows the 
analysis of these results in terms of corange and cotidal lines. 
The path of the M2 constituent tidal wave crest through the 
Gulf is described by the cotidal lines labeled in hours with 
zero being the time of the passage of the moon through the 
90th west meridian. 

The detail obtained in a 504 point grid considerably lessens 
subjectivity in these cotidal and corange analyses. Significant 
aspects of these analyses are: the well-developed amphidromy 
in the Northeastern corner of the Gulf, the presence of the 
greater amplitudes around the coast, the presence of a rather 
flat field of low amplitudes in the vicinity of the amphidromic 
point, and the relatively high amplitudes in areas of known semi- 
diurnal dominance. The path of the tidal wave crest through 
the Gulf is very well defined. 

Both the 16- and 504- point grid computer solutions 
produce corange and cotidal lines consistent with coastal data. 
Because one is based upon constant depth and the other upon 


actual depths the analyses differ considerably in some areas: 
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however, the regime shifts mentioned in Section 3 (say, 
semidiurnal to diurnal) are clearly indicated in both Figures 


4 and 5. 
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Yo Conclusions and Acknowledgements 

The three solutions obtained indicate that the numerical 
methods previously discussed are effective in predicting tides 
at sea, Based on the relatively quick convergence with actual 
depths of the 16 point system which has a grid distance of 
180 nautical miles, the convergence with constant depth of 
the 504 point system which has a grid distance of 30 miles, 
and the failure to converge with actual depths of the 50) 
point system, it is reasonable to assume that there exists an 
intermediate grid distance such that convergence may be ob- 
tained with maximum detail in the solution. Another approach 
is to obtain the maximum detail within the field by mathemati- 
cally smoothing the actual depth data to such a degree that the 
50% point system would converge. The danger here is that 
excessive smoothing could distort the data thus giving less 
realism in the solution. 

The cotidal and corange analyses of the constant and 
variable=depth models indicate that variations in depth are not 
solely responsible for the tidal regime shifts briefly described 
in Section 1; (see Figure 4), Additional studies may deliniate 
the relative importances of spacial forcing function variation, 
input current, basin shape, and basin depths upon tidal ampli- 
tude distributions and perhaps uncover other important variables. 

By removing certain topographical features of th Gulf floor 


from the data it would be possible to assess their effect upon 


=> 


tidal distributions. Similarly variations in input currents and 
the forcing function field could be examined for their effects. 
By adjusting Input currents only, under a condition of maximum 
realism and detail in all other data, to obtain the constituent 
distribution most compatible with coastal observations one would 
be able to evaluate the accuracy of present constituent tidal 
current data. 

The only geographical limitation in this technique is the 
size of the basin due to the beta-plane assumption. Thus all 
other seas comparable in size to the Gulf of Mexico are well 
suited for an investigation of this sort. 

This study and the writer have both benefitted greatly 
from the interest, guidance and considerable abilities of mem- 
bers of the faculty and staff of the United States Naval 
Postgraduate School. Principal among these are: Associate 
Professor T, Green, III for his encouragement, help, con- 
siderable effort and constructive guidance; Professor G. J. 
Haltiner for his counsel and personal example of scientific 
excellence; Associate Professor J. B. Wickham for his incisive 
and constructive criticism; R. D. Brunell for his personal 
interest, efforts in the programming and technical assistance; 
also Dr. G. B. Austin and Dr. G. Dowling of the U. S. Navy 
Mine Defense Laboratory, Panama City, Florida for their 
interest and assistance in selection of this topic; and Mrs. 
Julene L. Gainer for her patience, encouragement, understand- 


ing and recently acquired proficiency at the desk calculator. 
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Figure L. - RESULTS OF VARIABLE DEPTH 
COMPUTER SOLUTION: 


Corange lines (dashed) and cotidel lines (solid). 
Corange lines are labeled in centimeters. Cotidal lines 
are labeled in degrees denoting the phase of the tide 
relative to meridional passage of the moon through 
900W. 
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APPENDIX A 
SONI TE DIFPFERENCEFSEQUATIONS FOR INTERIOR 
AND ALL BOUNDARY POINTS ADAPTED FOR 
COMPUTER SOLUTION 
Equation (12) represents the finite difference form for 
an interior point. We have shown that at each boundary un- 
ique modifications of this equation are necessary to represent 
continuity and the boundary conditions. To avoid the necess- 
ity of separately establishing each equation within the program 
and to create a compact iteration scheme the following system 
was devised. 
Examination of all the equations led to a general form, 
(J-H)。 = а, ( SH) се Cb Som 
+Q2( 14) χ-ι + bz Pa 


+03 (T-H)a + O3( (2) vel 
+04 SH), Ба (Sad 


where the symbol 7 indicates that all terms on the right side 


are divided by C2 


E Aan eG (De +b ANR 
Ог 36909 - 61 (5-4) y. 


Е Оз (72) Yet 一 оз (FH) çe 
E G4 (2). e ba CXE- M), 
CU 
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These equations are set up to carry out Gauss-Seidel 
iteration, but modification for other solution methods could 


be made at this point. 


Now by indexing within the program the constants may 
be called to fit any of the previously derived equations satis- 
fying both dynamic considerations and boundary conditions. For 


interior points the constants are: 
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APPENDIX B 


DIFFERENCING EQUATIONS FOR VARIABLE 
GRID DISTANCE 


Given a relative grid arrangement of the sort, 


Y+\ 
A 


X+\ 





Ya 


it: is desirable to use only contiguous points in the differencing 
scheme because the use of other points would require a field 
of values and points outside the real grid. This in turn would 
negate the use of central differences in evaluating the second 
derivative. 


The basic forward differencing formulae were applied to 


о б) аз (Zr 2. fe e] 
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where Z is a dummy variable. This expression is quite involved 
and the stability of calculation about such a point is highly 


questionable. 
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APPENDIX D 


THE NUMERICABIR BESUAT ΟΕ πα ος ΕΙ 
CONSTANT BEPTFH БОЛО Ее 


POINT MAGNITUDE (cms) PHASE (Degrees) 


6. 1.3887000 —141.24518' 
ο 1.375051 -138.80641 
3. 1.3217907 -131.90854 
ly 1.553935; -120.89377 
Ee 1.2029531 -105.43202 
6. 1.2059750 -86 .843942 
Т 1.2334819 73. 334869 
8. 1.7517582 -172.5071 
9. 1.7194041 -172.30618 
πο, 1.7070080 -171.97713 
ei 1.6958950 -171.56044 
Nj 1.3027698 -110.26361 
13. 1.2262120 -137.87719 
E 1.1332772; -130.61788 
15. 1.0404123. -118.38878 
16. . 98849004 -99.810982 
I 1.0718649 -.75.709188 
15, 10555876 -53.916750 
19. 2 4490565 -39.475031: 
20. 4.1602358 -23.200350 
21% 7.4439120 -4.1130445 
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гез 
сз. 
el. 


Бо 


26. 
27. 
28. 
29. 
30. 
31. 
32. 
33. 
She 
う 5・ 
36. 
37. 
38. 
39. 
40. 
41. 
Ne. 
43 
hl. 
45. 
46. 


524827457 
1.6776942 
1.6361282 
ο ο... 
1 。5654 テ 88 
1 。556401 1 
1:5 
1. 5838993 
1.596388}, 
1.5907292 
125503097 
12405019722 
1 。25530 52 
1 。1279148 
1.0039504 
88153191 
79735342 
„8628115 
1.2173230 
1 。9446390 
2.9908028: 
1459681 
3. 7697796 
29949950 


1. 5900268 


P 


-2.1818875: 
-172.36824 
-172.15673. 
1171/6625. 
ー171 。39348 
-170.96102 
-170.38945 
-169.61342. 
-167.86595 
2165.26372: 
-161.59242: 
-156.27691. 
=18.15г22”7` 
ー142 。15830. 
-134.12128: 
-131.44061: 
-100.55519: 
-71.147404 
= ル 45 0294358 
-28 377334 
-15.251722 
-3 9740857 

1 . 7862483 

98396774 


-171.85187 


7. 
48. 
49. 
50. 
SE, 
52. 
53. 
She 
55. 
56. 
57. 
58, 
59; 
OD 
61. 
DES 
63. 
64. 
65. 
ος. 
C7. 
68. 
69. 
70. 
τα. 


1.526552 


OS ON 


1 , AMS2 


1.4258745: 
1 ROO 
1.424289}: 
1.4326960 
1.4367689 
1.4264928: 
1.3910 
1.315395% 
1. 1943604. 
1.07521430 
. 94189596 
79313253 
・64595142 
・59415964 
ӘСІ 
Lao 
1.8856044 
2. uu22500 
2.5209018 
2.53860145 
1.583 56 
1.5106958 
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-171.55308 
-171.29280 
-171.10465. 
-170 . 47495 
-169.98678. 
-169.33126 
-168.36735 
-166.84673. 
-164.70118 
-161.82801 
-158.02929 
-153.24541 
-148.40719 
=1}2.03557 
-132.02672 
“11855350 
ΙΤ 
-44.336330 
-22.528501 
-9.01,38157 

1.3021725 

10.436640 
20.977383 
210135588 
SUE а 


72o 
73. 
Th o 
75° 
76. 
TTo 
78. 
79: 
80. 
81. 
82. 
83. 
δὴ ο 
85. 
86. 
87. 
88. 
89. 
90. 
91 。 
92. 
93. 
Qh o 
95・ 
96. 


14459243 
1.3873289 
13523482 
1.3203357 
1.301490) 
1. 3007594 
1.3081634 

1.3162539 
1.3232527 
1.3207365 
1.3008327 
55 
1.179141) 
1.0871412 

968,441 

. 82084260 

-635086931 

42310620 

32467912 

. 55898858 

.93149353 
1.295355 
1.99208 
1.8510571 


1.4426962 
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~170 .90635 
-170 .62673 
-170.37487 
-170.0017), 
-169.57599 
-169.12103 
-168.52707 
-167.68322 
-166.52922 
-165.02452 
-163.15918 
-160.91726 
-158.41437 
—155. 84092 


-152.64950 
-148.05767 


-139.73622 
-119.25154 
-66 . 502209 
=18 . 151299 
2 .2351259 
15.742440 
29.348534 
43.251027 


-170.27041 


97. 
98. 
99. 


100. 


1.3801700 
1.31Ə55B4. 
1 。2602636 
1.2306716 
1.2046669 
1.1925103 
1.1990822 
1.2142370 
1.220782), 
1.2495502 
1.2644 726 
1.2694225 
1.258582 
1.2247704 
1.1'738689 
1.0947432 
.98741067 
‚8412120 
・64327279 
・40957519 
‚16316727 
ο 68Ε4. 
«60679389 
. 98720807 


1 。 テ 865882 
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-170.11656 
-169.87776 
-169.60957 
-169.39582 
-169.06696 
-168.71396 
-168.38716 
-167.98880 
-167.44851 
-166.79710 
-166.04208 
-165.22816 
-164.40579 
-163.71041 
-163.24897 
-162 .92580 
-162.80156 
-162.89192 
-163.50925 
-167.70662 

160 。04043 


75.916858 
6745911211 
74.397350 
81.249956 





122. 
123- 
124. 
1253 
126. 
127- 
128. 
129. 
130. 
131. 
132. 
133. 
134 . 
195% 
136 。 
137 。 
138. 
139. 
140. 
ц. 
Wwe. 
143. 
Wd 
145. 
146. 


1. 3192293 
1.2570372 
1.1937553 
1:14,50 593 
1 .1211817 
1.1020564 
1. 0981462 
1.145761 
1.15354 
1.1725985 
1. 2102442 
1.2500523 
1.2873929 
129199386 
1.3320425 
1 。3344956 
1.3135021 
1.2716266 
1.2020871 
1.0981960 
29868950 
. 88220144 
. 8666759 
。98653169 
1。2799964 
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-168.98529 
~168 .94400 
-168. 734,72 
=168 。49526 
=168。34803 
=168 。10229 
-l67 . 881417 
=167 。77086 
=167 。67201 
-167 . 54610 
-167 .46968 
=167 。47233 
=167 。60798 
=167 。92557 
-168. 5213} 
-169 .1,7107 
-170 . 771,50 
172 . 57h 
-175 0535] 
-178.68367 
175.47563 
165 .61313 
150 .92426 
136.05371 
196 «90135 


1,7. 


148. 


149. 
150. 
151 。 
152. 
153. 
154 
155. 
156. 
157. 
158. 
159. 
160. 
161. 
162. 
163. 
164 . 
165. 
166. 
167. 
168. 
169. 
170. 


171. 


1. 7474118 
12092843 
1.161983 
1,0849542 
1。0403817 
1.02202 
1.0109213 
1.0166113 
1。0451478 
1.0872336 
1.1373998 
1.1995002 
NANE 
1.370828 
] . 4228629 
1.49577548 
1. 5614789 


7 1。6100079 


1. 0439044 
1.656183] 
1.610146 

1.6223832 
1.6032476 
1.6392330 
l. 7547253 
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121. 6906} 
-167 。64318 
-167.67550 
-167.48544 
-167.27501 
-167.21450 
-167.08640 
-167.05710 
-167.22189 
-167.49308 
-167 . 83048 
-168 。 32647 
-169.00223 
-169.88820 
-171.00051 
-172.38357 
-174.05059 
-175.98955 
-178.23670 
179.1847 
176 .18332 
172.66604 
168.28565 
161.61001 


155.14506 


172% 
173. 
174. 
175 。 
176 。 
177. 
178. 
179. 
180. 
181. 
182. 
183. 
18). 
185. 
186. 
187. 
188. 
189. 
190. 
191. 
192. 
193. 
194 ο 
195. 
196. 


20045501 
2.3098791 
3.05812 
1.10866753 
1.044;2983 
-98475467 
-94421575 
093236154 
92907261 
o 94545533 
.987'70361 
1.0474125 
1,1201448 
1. 211263 
123187259 
lo 4407064 
la 5713587 
1. 7083588 
1. 8462118 
1. 9728861 
2.0883663 
261825526 
2.29039 
2.2967387 
233770167 
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150 .74902 

149 o 75370 

155.39824 
=166 。33404 
=166 。35663 
=166 。1313 ル 
=165 。92555 
=165 。95960 


=165 。97488 


-166. 17686 
=166 „67282 
=167 。33990 
=168 。13757 
169 .144 57 
=170 。35560 
=171。76804 
=173 。36616 
=175.15459 
„177 11356 
-179 . 21124 
178.57807 
176.30378 
17l. 02875 
172 .02866 


171, 09147 


Torm 
198. 
199. 
200. 
201. 
202. 
203. 
204. 
205. 
206. 
207. 
208. 
209. 
OE 
eid 
eld 
ei 
214. 
215. 
ZUR 
217% 
eU. 
219. 
220. 


フン 


2.4036320 
2. 5284900 
2. 7636824 
3.1167505 
3. 5327726 
1.0129899 
. 94788364 
0 89049944 
‚85424915 
. 84869031 
85389577 
68147957 
«93842593 
1.0173279 
1.1149818 
1.2381205 
1.3866221 
1.5595657 
1.7520227 
1 .961617 ル 4 
2.1814858 
23965146 
2.6017553 
2.7802197 


2.9135588 


δὴ 


165481088 
162.15293 
159.89052 
158. 85404 
157.76749 
-165.13166 


-164.99387 


-164.63719 
、 =164。38986 


-164.51702 
-164. 69366 
-165.15931 
-166.01340 
-167.08388 
-168.30477 
-169.73308 
-171.33510 


-174.95051 
-176.92430 
-178.97981 
178.91438 
176779859 
Nene 
17274419 


ERËS 
гави 
gel 
2250 
220% 
RRE 
228% 
229. 
230. 
2347 
2381! 
295% 
gal. 
2 55 
236. 
237. 
ве» 
239. 
ΞΘ. 
eM. 
— 
AUS 
eh. 
ぐう < 
26. 


3.0161880 
3.0826239 
351604651. 
3.2728319 
3.189903), 
3.8443713 
1, . 3693820 
5825 921 
1014823815 
.91611772 
. 85362642 
3709/6993 
«76826352 
. 76900960 
8258238 
„82125881 
.89298309 
99155335 
19252134 
Wal 
1 -4649441 
1 。6935393 
1 。95 ル 46477 
2。2r64576 


る 。5630418 


БЕ 


170 954854 
169。29360 
167.06122 
165.33070 
164. 56022 
164.85080 
166.24802 
169.77120 
-165.60600 
-164.13332 
-163. 50215. 
=162 。89356 
-162.55455 
10277507 
~163.13041 
-163.88891 
-165.12379 
-166 . 58995 
-168.18580 
=169 。94772 
-171.82287 
-173.77752 
-175.78511 
-177.84262 


179. 99587 


BUT « 
248. 
249. 
250. 
ο 
252% 
253; 
25%. 
255. 
256. 
257. 
258. 
259. 
260. 
261. 
262. 
263. 
261. 
265. 
266. 
267. 
268. 
269. 
270. 


i 


28795196 
3 。 1896662 
3.4633445 
3.6690512 
3.8103315 
3.8769698 
3.9217889 
3.9790625 
4, 1409762 
4 «4609625 
469953985 
5. 8880951 
2656523 
8.1811000 
„90562088 
‚81687892 
. 76041280 
. 71118680 
・68455013 
69116844 
„71238788 
‚76123619 
・84674989 
. 96413752 
1.111 


66 


177.95661 
17525195 
173.86999 
171.98936 
178 731615 
168.82913 
167.55483 
166. 78406 
166. 88418 
167.91345 
169.71161 
171.96057 
17321906 
17554859 
-163.4260% 
-162.22822 
-161.)8515 
-160.67153 
-160.23011 
-160 . 56392 
-161.12161 
-162.21030 
-163.85739 
-165.71847 
-167.65178 


272. 
273. 
27k, 
215 
276. 
277. 
276% 
279. 
280. 
281. 
282. 
283. 
284 . 
285. 
286. 
287. 
288. 
289. 
290. 
291. 
292. 
293. 
29%. 
295. 
296. 


1 。3048910 
1.51,28370 
1.8304559 
2.1665757 
2.551392 
2.9773257 
3.4219972 
3.8637190 
le 257099 
4. 5479984 
47205917 
e 7548153 
4. 7096263 
Ņ -6451493 
4 。6905389 
Іі. 9106693 
5。3524610 
6.132279) 
7.3026670 


8. 8962829 


¿80854982 
. 12256166 
。66983155 
。62448956 
«60207424 


67 


-169.68291 
-171.75315 
-173.83722 
-175.92528 
-178.03303 
179.83285 
177.67307 
175.52600 
173.15832 
171.53646 
169.86121 
168.47200 
167. 52009 
167.17105 
167.78207 
169.35573 
171.62019 
17) .05608 
176.00718 
177.29166 
-161.11039 
-159.74800 
-158.829),5 
=157 。74859 
=157 。16038 


297. 
298. 
299. 
300. 
301. 
308% 
303. 
304. 
3958 
306. 
307. 
308. 
309, 
310. 
Зі 
3154 
313. 
31h. 
3158 
316% 
317. 
318, 
319. 
320. 


3218 


・61349231 
・64091 ル 441 
. 69805652 
79510355 
.92885725 
1.101186) 
143251186 
1.6076159 
1.9553803 
μη 
2.8598067 
3.4176876 
l, «0214149 
4 。6391439 
5.1997527 
5。60825 ル 4 ル 
っ ・8143937 
5e 7769834 
595581712 
5.2677903 
51227911 
5.1654102 
5.458114,57 
6 。0624303 
6.981911) 


68 


-157.63016 
~158.42782 
-159.90782 
-162.02933 
-161,.31382 
-166.57697 
-168 . 84837 
-171.07635 
-173.25640 
-175.40321 
-177 . 55745 
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